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O ! ABSTRACT 

We study the effects of mergers on the structural properties of disc-like systems by 

, using Smooth Particle Hydrodynamical (SPH) numerical simulations in hierarchical 

^ ' clustering scenarios. In order to assess the effects of mergers on the mass distributions 

00 I we performed a bulge-disc decomposition of the projected surface density of the sys- 

i tems at different stages of the merger process. We assumed an exponential law for the 

' disc component and the Sersic law for the bulges. We found that simulated objects at 

I z = have bulge profiles with shape parameters n ~ 1, consistent with observational 

^"-^ i results of spiral galaxies. The complete sample of simulated objects at z = and z > 

' shows that n takes values in the range n ~ 0.4 — 4. We found that secular evolution 

, tends to produce exponential bulge profiles, while the fusion of baryonic cores tends 

• to increase the n value and helps to generate the correlation between B/D and n. We 

Qh| found no dependence on the relative mass of the colliding objects. Our results suggest 

^ . that mergers, through secular evolution and fusions, could produce the transformation 

' of galactic objects along the Hubble sequence by driving a morphological loop that 

I might also depend on the properties of the central galactic potential wells, which are 

Ci • also affected by mergers. 
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1 INTRODUCTION 

The origin of the Hubble sequence is still a controversial is- 
sue. In particular, spiral galaxies seem to evolve along this 
phenomenological classification although the physical mech- 
anisms behind these morphological changes are not fully un- 
derstood. The structural parameters, the gas abundances 
and the star formation activity vary along galaxies in the 
Hubble sequence in the sense that disc-dominated systems 
are also the more gaseous ones, experiencing on-going star 
formation. The large database of galaxies gathered in the 
last years have allowed to get more detailed information on 
the properties of different morphological type galaxies in the 
local Universe. 

Observational results show that when a double expo- 
nential decomposition is applied to the sufarce luminosity 
density of late-type spirals, the scalelengths of the bulge 
and disc components seem to be restricted to a certain 
value < rb/rd >~ 0.10 (Courteau, de Jong & Broeils 1996, 
hereafter CdJB96). This restriction in the range of possi- 



ble scalelengths for spirals has been interpreted as a proof 
for secular evolution to be the responsible mechanism for 
bulge formation from an already in place disc structure. Re- 
cently, MacArthur, Courteau & Holtzman (2002, hereafter 
MCH02) studied a larger sample of late-type spirals and 
carried out a bulge-disc decomposition assuming a Sersic 
law for the bulge surface luminosity density. These authors 
found a continuous distribution of the shape parameter of 
the bulge (n « 0.2 — 2) with a maximum at n ~ 0.90. They 
also found a restricted range of values for the scalelengths 
of the bulge and disc when the shape parameter is n « 1, 
< rh/rd >~ 0.13±0.06 in the ii-band. Previous works have 
also found that bulges of spirals can be fitted by a Sersic law 
with different shape parameters (e.g., Andredakis, Peletier 
& Ballcels 1995, hereafter APB95; Khosroshahi, Wadadekar 
& Kembhavi 2000, hereafter KWKOO; Graham & de Block 
2001, hereafter GdBOl). APB95 and GdBOl also found a 
correlation between the luminosity bulge-to-disc ratio B/D 
and the shape parameter of the bulge, which also correlates 
with morphological type. These observational results suggest 
a possible connection between the formation of the bulge and 
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disc components. They also support the idea that n could be 
used as a good indicator of position in the Hubble sequence. 

Several theories have been developed to explain the for- 
mation of the disc and bulge components. In the case of the 
disc systems, the standard model is based on throe hypothe- 
sis: the angular momentum is acquired through cosmological 
torques (Peebles 1969), baryons and dark matter have the 
same specific angular momentum content (J), and this spe- 
cific angular momentum is conserved during the collapse and 
coohng of baryons (Fall & Efstathiou 1980, hereafter FE80). 
This model has been successful in reproducing several ob- 
servational results in analytical and semi-analytical mod- 
els (e.g., Dalcanton, Spcrgel & Summers 1997; Mo, Mao & 
White 1998). However, serious problems arose in numerical 
simulations where an important angular momentum trans- 
fer from baryons to the dark matter haloes during merg- 
ers was detected, breaking the condition of J conservation. 
Dommguoz-Tcnroiro, Tissera & Saiz (1998, hereafter DT98) 
showed that a disc-like structure with observational coun- 
terpart (Saiz et al. 2001, hereafter SOI) can be built up 
if a compact stellar bulge is allowed to form without de- 
pleting the gas reservoir of the system. These stellar bulges 
provide stability to the gaseous disc systems which are ca- 
pable of conserving a non-negligible fraction of its angular 
momentum during violent events. Supernova energy feed- 
back could also contribute to the formation of the disc com- 
ponent, regulating the star formation rate and preventing 
early catastrophic depletion of the gas into stars. Probably 
both mechanisms, the formation of a compact stellar bulge 
which assures the axisymmetrical character of the potential 
well and energy feedback, work together in nature to allow 
the formation of spiral galaxies (e.g. Weil, Eke & Efstathiou 
1998). 

The formation of bulges is a more complex task since 
several mechanisms could be acting together such as mono- 
lithic collapse (Gilmore & Wyse 1998), mergers (Kauffman, 
Guiderdoni & White 1994) and secular evolution (Pfenniger 
& Norman 1990). In general, analytical models and pre- 
prepared simulations have focused in one or two of them 
at the time (e.g., van den Bosch 1999; Aguerri, Balcells & 
Pclcticr 2001, hereafter AOl). These models have been suc- 
cessful in explaining some properties of spiral galaxies, de- 
spite their approximations. 

The current paradigm for the formation of the struc- 
ture favors a hierarchical clustering scenario where the struc- 
ture forms by aggregation of substructure. Hence, a galactic 
object experiences the effects of collapse, merger, interac- 
tion and probably secular evolution in a non simple fashion, 
which can also strongly depend on rcdshift. In particular, 
violent events (i.e. mergers, interactions) can have impor- 
tant effects on the internal properties of the objects, such as 
their mass distribution (Mihos & Hernquist 1996) and star 
formation activity (Tissera et al. 2002, hereafter T02). 

Minor and major mergers of disc systems with satellites 
have been extensively studied. However, most of these works 
do not consider the presence of the stellar bulges. The rel- 
evance of this component in the overall evolution of a disc 
system during a merger has been firstly pointed out by Mi- 
hos & Hernquist (1994) in a study of pre-prepared mergers. 
Recently, AOl studied the effects of mergers on the struc- 
tural parameters of bulges of disc-like systems by assuming 
that initially bulges form with an exponential profile. These 



authors found that coUisionless mergers produce a migration 
to higher n and that this effect is proportional to the relative 
mass of the colliding systems. The simulations used in this 
analysis were non-cosmological pre-prepared ones where on 
going star formation and gas dynamics were not included. 
Recently, T02 show for the first time in cosmological SPH 
simulations how the central mass concentration in disc-like 
systems can grow by collapse, merger and seculax evolution. 

These numerical results suggest that a key point in the 
formation of both bulges and discs in hierarchical clustering 
scenarios seems to be mergers. Mergers have been tradition- 
ally suggested as a possible mechanism to drive the Hubble 
sequence since they can modify the dynamical and astro- 
physical properties of galactic objects. Hence, the question 
would be how, in hierarchical scenarios, mergers work to 
shape disc-like objects and if the properties of these objects 
resemble those of current spirals. In this paper we will focus 
on the analysis of the mass distributions of disc-like systems 
and how they are modified during mergers, paying special 
attention to the comparative study of the effects of secular 
evolution triggered by tidal fields and of the actual colli- 
sion of the baryonic clumps. In order to assess the effects 
of mergers on the mass distributions we perform a bulge- 
disc decomposition of the projected surface density of the 
systems at different stages of the merger process. We then 
study how such different structural parameters, including 
the shape parameter defined by Sersic (1968), evolve during 
the orbital decay phase and fusion of the satellite. 

Our simulations include the effects of gravitation, hy- 
drodynamics, cooling and stax formation in a cosmological 
framework. Hence, the set of mergers that we analyse axe 
given by the particular evolutionary history of each galactic 
object. This is a crucial point since the merger parameters 
and physical characteristics of the colliding objects are not 
ad-hoc choices but result from the consistent formation of 
the structure in a hierarchical scenario. The drawback of 
our approach is a lower numerical resolution compared to 
those used in studies of pre-prepared mergers. Assessment 
of possible numerical problems are discussed throughout the 
paper. 

In section 2 we present the analysis of the simulations. 
In section 3 we discuss the results. Section 4 summarizes the 
conclusions. 



2 MODELS AND ANALYSIS 

We run cosmological SPH simulations which include hydro- 
dynamics and star formation (Tissera, Lambas & Abadi 
1997). The simulated boxes represent typical regions of a 
standard Cold Dark Matter (CDM) universe of 5 Ir^ Mpc 
side and 64^ particles {Q = 1, A = 0, h = 0.5). We 
assume a baryonic density parameter of fib = 0.10. All 
baryonic and dark matter particles have the same mass, 
Mp — 2.6 X 10*Mq . We have used a gravitational softening of 
3.0 kpc and a minimum hydrodynamical smoothing length 
of 1.5 kpc. We performed three simulations with different 
realizations of the power spectrum (cluster normalized) and 
the same cosmological and star formation parameters. The 
simulated volume is the result of a compromise between the 
need to have a well-represented galaxy sample and enough 
numerical resolution to study the astrophysical properties of 
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the simulated galaxies. We are confident that since we focus 
our analysis on small scale processes such as mergers and 
interactions, scale fluctuations of the order of the box size 
will have no significant effect on such local processes. Note 
also that we are interested in the effects that mergers have 
on the mass distributions of the galactic objects analysing 
them as individual events and not in connection with their 
evolution or environment. 

The star formation algorithm used in these models is 
based on the Schmidt law (see Tissera et al. (1997) for de- 
tails). Cold gas particles are eligible to form stars if they 
are denser than a certain critical value and satisfy the Jeans 
instability criterium (Navarro & White 1994). Gas particles 
are checked to satisfy these conditions at all time-steps of 
integration. Hence, as the gas cools down and is gathered at 
the centre of dark matter haloes, it is gradually transformed 
into stars according to the particular history of evolution of 
each galactic object. Only one free parameter has to be fixed: 
the star formation efficiency, which links the dynamical time 
of the gas cloud represented by a particle with its star for- 
mation timescale. We have used the value adopted by SOI 
since it is adequate to reproduce disc-like structures with 
observational counterparts at z = 0. These authors used a 
low star formation efficiency which allows the formation of a 
stellar bulge that assures the axisymmetrical character of the 
potential well but without exhausting the gas reservoir (see 
also DT98). As a consequence, disc-like structures can be 
formed, although they remain mainly gaseous. Conversely, 
observed discs are mostly stellar, but they inherit the struc- 
tural and dynamical properties of the gaseous discs out of 
which disc stars are formed. Hence, for the sake of compari- 
son, the values of the parameters describing those properties 
can be reliably estimated from the simulated gaseous discs. 
In this paper, we will determine and use these parameters 
to study the effects of mergers on the mass distribution of 
the simulated galactic objects. 

At 2 = 0, we identify galactic objects at their virial ra- 
dius, analysing only those with more than 4000 dark matter 
particles within their virial radius in order to diminish nu- 
merical resolution problems. Within each galactic object, a 
main baryonic clump is individualized which will be, here- 
after, called the galaxy-like object (GLO). 

The selected GLOs have very well-resolved dark matter 
haloes which provide adequate potential wells for baryons to 
collapse in. This fact assures a reliable description of the gas 
density profiles (see Steinmetz & White 1997), which allows 
us to follow the star formation history of the GLOs (see 
Tissera 2000). With this strong restriction on the minimum 
number of particles, the final GLO sample at 2 = is made 
up of 12 GLOs with virial velocities in the range 140-180 
km-s^^. 

We follow the evolution of the selected GLOs with look- 
back time {t{z) = 1 — (1 -|- z)~^^^) identifying mergers and 
starbursts (SBs). We then construct their mergers trees and 
star formation rate (SFR) histories (because of the restric- 
tive SF efficiency used, these SFR histories are mainly those 
of the bulge components). During a merger event, the pro- 
genitor object is chosen as the more massive baryonic clump 
within this merger tree, while the minor colliding baryonic 
clump is referred to as the satellite. Following T02 we study 
only those mergers that are directly linked to stabursts in 
the SF history of each GLO. These authors found that dur- 



ing the orbital decay phase of some merger events, early 
gas inflows which trigger star formation can develop. In this 
case, if the gas reservoir of the GLOs is not exhausted dur- 
ing these starbursts, second bursts are induced at the ac- 
tual physical contact of the baryonic clumps. The induction 
of early gas inflows could be directly linked to the proper- 
ties of the total potential wells of the systems in agreement 
with previous results reported by Barnes & Hernquist (1996) 
and Mihos & Hernquist (1996) in pre-prepared simulations. 
When no gas inflows are triggered during the orbital decay 
phase (ODP), only one SB is detected when the two bary- 
onic clumps collide. These events will be called single SBs 
(SSBs), while those where two SBs are detected will be called 
double SBs (DSBs). We studied a total of 18 merger events 
with merger parameters settled by the cosmological model. 
Among them, 11 are classified as DSBs and 7 as SSBs. Fig. ^ 
shows an example of both types of SBs. We have also plotted 
the distance between the centres of mass of the progenitor 
and the satellite since the redshift they share the same dark 
matter halo. The merger event will be determined by four 
redshifts of reference, za'- the beginning of the first bursts 
in DSBs and the ODP in SSBs, zb'- the end of both the first 
bursts in DSBs and the ODP in SSBs, zc'- the beginning 
of the second SBs in DSBs and the only bursts in SSBs (in 
this case zb = zc), and zu'- the end of the second bursts 
in DSBs and of the only one in SSBs. These redshifts have 
been depicted in Fig. ^. Note that the ODP, defined as the 
time since two objects share the same dark matter halo to 
the fusion of their main baryonic cores, is determined by za 
and Zc ■ The zb redshift is used in the case of DSBs since it 
establishes the end of the first SBs. 

In order to quantify the changes in the baryonic mass 
distributions, we carried out a bulge-disc decomposition of 
the projected mass surface density of the GLOs at z = 
and of the progenitors at each of the redshifts that define 
a merger event. We assumed an exponential profile for the 
disc component and the Sersic law for the bulge central mass 
concentration: 

Ed = E5exp(-(r/rd)) (1) 

Eb = ESexp(-(r/rb)^/") (2) 

where and r-b are the disc and bulge scalelengths, E^ and 
Eb the corresponding central mass surface densities and n 
the bulge shape parameter. 

For each GLO at z = and its progenitor objects, we es- 
timated the total angular momentum (j) and projected the 
mass distribution onto the perpendicular plane defined by j. 
We then integrated the projected baryonic mass in concen- 
tric cylinders of radius r since, as it has been proved by SOI, 
it is numerically more convenient to fit the integrated mass 
than the surface density. In order to diminish the effects of 
numerical resolution, the fits were made from r = 1.5 kpc 
to r « 30 kpc, except in those cases where the mass distri- 
butions of the progenitors were strongly perturbed by the 
entrance of the satellites. 

The total integrated baryonic mass of the GLOs 
(AfcLo) is estimated, on one hand, as the baryonic mass 
within r — 30 kpc, and will be associated with the luminous 
galaxy. On the other hand, Mglo is the result of the disc 
(Md) and bulge (A'lt) mass contributions: 

AfGLO = A/d + Mb (3) 
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Figure 1. Star formation rate (SFR) for typical galaxy-like ob- 
jects during a double starburst (a) and a single one (b) as a func- 
tion of look-back time (t(z)). The solid lines represent the am- 
bient star formation components and the dashed lines depict the 
distance between the mass centers of the progenitors and satel- 
lites (D). The four redshifts of interest during a merger event have 
been also plotted. 



Mb = 2Tv'EZrlnr{2n) 



(4) 
(5) 



where F is the complete gamma function. Since Mglo is 
known, there are four free parameters to be fitted: rd, rb, 
Eb (or Ed) and n. From these fitted parameters we analyse 
possible correlations with the dynamical and astrophysical 
properties of the GLOs and their dark matter haloes. In par- 
ticular, we will study how the structural parameters change 
during merger events. 

Note that we work with masses instead of luminosities, 
hence, mass-to-light ratios for the bulge and disc compo- 
nents will have to be adopted in order to be able to com- 
pare our results with observations. To this end we use dif- 
ferent observational results from APB95 (30 SO-Sbc galaxies 
in the A'-band), CdJB96 (173 spiral Sb-Sc galajcies in deep 
J?-band), Moriondo, Giovanelli & Haynes (1999, hereafter 
MGH99, 25 Sa-Scd spiral galaxies in the A'-band), KWKOO 
(27 early-type disc galaxies in the A'-band) and MCH02 (47 
late-type spiral galaxies in the i?-band). APB95, KWKOO 
and MCH02 used bulge-disc decompositions of the luminos- 
ity surface density profiles, allowing the bulge shape parame- 
ter to vary. CdJB96 used double exponential fittings (n = 1). 
Finally, MGH99 made fittings with different n values find- 
ing that most of their bulges were better described by n = 1 
(except for two with n = 2) . Note that all these observations 
correspond to galaxies at z « 0. 
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Figure 2. Disc scalclcngth (rd) as a function of bulge scalelength 
(rb) for simulated objects at z = (filled pentagons). Labels with 
the shape parameters (n) of the simulated bulges are shown. Ob- 
servational data from CdJB96 (crosses), MGH99 (n = 1 crossed 
squares, n = 2: open squares), KWKOO (n < 1.5: dashed tri- 
angles, n > 1.5: open triangles) and MCH02 (n < 1.5: crossed 
circles, n > 1.5: open circles) are also included. 



3 RESULTS 

3.1 GLOs at z = 

As we mentioned before, SOI showed that disc-like objects 
formed in numerical simulations where stellar bulges were 
allowed to form have structural parameters and specific an- 
gular momentum content similar to those of observed spi- 
ral galaxies. For the sake of completeness and to show that 
these findings are valid for our GLOs, we will first discuss 
their properties at 2: = 0. The structural parameters and 
the circular velocity {V'^ = GM(r)/r) a.t r — 2.2rd allow 
the classification of our objects as intermediate spirals (100 
km-s"^ < V2.2 < 180 km-s~^; 4 GLOs), bright spirals (V2.2 > 
180 km-s^^ and rd > 5.25 kpc; 2 GLOs) and compact bright 
spirals (V2.2 > 180 km-s"^ and rd < 5.25 kpc; 6 GLOs). This 
is the result of our strong condition on the mass of GLOs in 
order to work with those better resolved. 

In Fig. ^ we show the relation between the disc and 
bulge scalelengths. Each GLO has a label that specifies its 
n shape parameter. We include the observational data from 
CdJB96, MGH99, KWKOO and MCH02. Note that for in- 
creasing n parameters, the scalelengths tend to move to the 
upper envelope of the distribution. This effect is also seen 
in the observed scalelength distribution. We have adopted 
ri = 1.5 as a general cut-off value to segregate the samples 
between those with exponential and non-exponential profiles 
for the bulge components. 

In Fig. ^ we show the normalized histrogram of n for 
GLOs at 2; = (a), and the corresponding distribution for 
the observations of APB95, KWKOO and MCH02 (b), nor- 
malized to the total number of galaxies in each sample. 
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Figure 3. Histogram of the sliapc parameter of the bulges for simulated objects at 2 = (a) and observations from APB95 (dashed 
line), KWKOO (dotted line) and MCH02 (solid line), normalized to the total number of galaxies in each sample. 



The combined observational sample covers a morphologi- 
cal range from Sa to Sc galaxies, but this is not a consis- 
tent sample. We show it in order to broadly compare the 
range of n values obtained from the simulations with obser- 
vations. The average n values over the observed subsamples 
are < uapbss >= 2.97±1.29, < nKWKOo >= 2.88±0.94 and 

< »iMCH02 >= 1.02 ± 0.43 while the simulated sample has 

< Usim >= 1.23 ±0.40. The simulated distribution is consis- 
tent with the observations of MCH02 at one cr-level, which 
suggests that our objects are more similar to late- type galax- 
ies. However, < nsim > value also agrees at 3 cr-level with 
the corresponding mean values of the early-type samples. 

A key point in the formation of discs is the conservation 
of the specific angular momentum (FE80). Our SF algorithm 
has been implemented in such a way that it has allowed the 
formation of stellar bulges but without exhausting the gas 
reservoir (DT98). Hence, gaseous discs have been able to 
regenerate after mergers. In order to compare the angular 
momentum (j) of the three mass components: bulge, disc 
and dark matter halo, with observations, we plot in Fig. 
^ their specific angular momentum (J=j/M) content as a 
function of mass (A/). The mass and angular momentum 
of bulges have been estimated with stars within r < rcff , 
while those of the gaseous discs have been calculated within 
r < 3.2 rd. These criteria have been adopted following those 
used by observers. The mass and angular momentum of the 
dark matter haloes have been estimated at the virial radius. 
In this figure, we have also included two boxes that depict 
the observational region covered by spirals and ellipticals as 
given by Fall (1983). From this plot we can see that the 
disc components and the dark matter haloes have compara- 
ble specific angular momentum contents as predicted by the 
standard disc formation model of FE80. Conversely, stellar 
bulges have been formed from material that has lost most 
of its angular momentum. 

From these results we conclude that the simulated 
GLOs have structural and dynamical parameters that sta- 
tistically resemble those of current normal spirals at z = 0. 
In the following section we will study the role played by 
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Figure 4. Specific angular momentum (J) for the gaseous discs 
at 3.2r;j (filled squares), the stellar bulges at r^ft (open stars) and 
dark matter haloes at the virial radius (filled triangles), as a func- 
tion of mass (Af ). The solid (dotted) box depicts the observational 
area for spirals (ellipticals) from Fall (1983). 



mergers in the determination of these properties and try to 
understand the origin of correlations among the structural 
parameters such as B / D versus n. 



3.2 Analysis of progenitors during merger events 

We carried out the bulge-disc (B-l-D) decomposition of the 
progenitor objects of the GLOs analised in the previous sec- 
tion during merger events at the four redshifts of interest: 
za, zb, zc and zq. In this section we will study how their 
structural parameters change during these processes. 
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Table 1. Mean values of bulge scalelength (rb), disc scalelength 
(rj) and bulge shape parameter (< n >) for single (S) and double 
(D) SBs for the redshifts of interest. 







rb [kpc] 


rd [kpc] 


< n > 


2A 


s 


0.46±0.17 


6.2±0.7 


1.5±0.2 




D 


0.31±0.11 


4.9±0.6 


1.8±0.2 


ZB 


s 


0.57±0.16 


5.6±0.7 


1.3±0.2 




D 


0.37±0.16 


5.9±0.4 


1.6±0.2 


ZC 


S 


0.57±0.16 


5.6±0.7 


1.3±0.2 




D 


0.53±0.12 


5.0±0.4 


1.2±0.1 


2D 


S 


0.52±0.12 


5.0±0.7 


1.4±0.3 




D 


0.37±0.08 


5.4±0.6 


1.4±0.2 



Note; The errors correspond to the bootstrap method, a"bt- 
3.2.1 The shape parameter 

In Fig. ^(a) we show the rd versus rb relation for all pro- 
genitors during merger events and at z = 0. We have also 
included the observational data from a B+D decomposi- 
tion from MGH99, KWKOO and MCH02, and from a double 
exponential decomposition by CdJB96. The combined ob- 
served sample of spirals shows a certain correlation between 
rd and rb, although some of them have, for a given disc 
scalelength, a smaller bulge one than that expected from 
this relation. We also note that the scalelengths of our sim- 
ulated disc structures show a similar behaviour. In order to 
individualize which observed and simulated objects belong 
to each of these two distributions of scalelengths, we seg- 
regate them according to their shape parameters as shown 
from Fig. |^(b) to Fig. |^(d). It is clear that as galaxies and 
simulated objects with large shape parameters are taken out, 
the relation for both, observed and simulated scalelengths, 
gets better defined. 

We also note that the rb versus rd relation for GLOs and 
galaxies with bulges with n > 1.5 shows no trend at all. It 
is encouraging that the observed and simulated disc struc- 
tures show scalelength distributions with similar patterns, 
for both n ~ 1 or n 7^ 1 bulge profiles. In the case of expo- 
nential bulges (0.9 < n < 1.1), the correlation is well defined 
with < rb/rd >= 0.14 ± 0.05 and < n^/n 0.13 ± 0.06 
for observations and simulations, respectively. Table ^ sum- 
marizes the mean values of rb, rd and n for GLOs at the 
four z of interest. 

Similary to Fig. ^(a) we display in Fig. |^ the histogram 
of shape parameters for all progenitors at all studied red- 
shifts. Hence this distribution gathers the information on 
the n parameters of different progenitors at different stages 
of evolution. From the comparison with that shown in Fig. ^ 
corresponding to the bulges of the final structures at z — 0, 
we deduce that in the past of current GLOs, their progen- 
itors went through stages where their bulges had n > 1.5 
parameters (~ 35 per cent of the total sample). Meanwhile, 
GLOs at z = show a trend to have n ~ 1, with less than 
10 per cent being described by n > 1.5. 

In order to assess the effects of mergers on the shape 
parameters of the objects and to look for possible differ- 
ences between the effects of secular evolution and fusions, 
we calculate the percentages corresponding to the increase 
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n 

Figure 6. Total distribution of bulge shape parameters for the 
progenitors during merger events. 



Table 2. Percentage of objects which increase (t n)/decrease 
(J, n) their shape parameter during the ODP {z\-zc), the fusion 
of the baryonic clumps (zq-zy)) and the whole merger process 
{zq-zd), for single (S) and double (D) starbursts. 







ZA-zc [%] 


zc-zn [%] 


ZA-ZD [%] 


s 


Tn 


42.9 


57.1 


57.1 




i n 


57.1 


42.9 


42.9 


D 


Tn 


18.2 


63.6 


45.5 




In 


81.8 


36.4 


54.5 



and decrease of the shape parameter during the intervals 
za-zc, zc-zu and za-zo for SSBs and DSBs. The results are 
shown in Table |[ Note that we are taking the interval za- 
zc as the total period corresponding to the ODP. For both 
single and double SBs, we find that bulges tend to decrease 
their n parameters during the ODP, and that this effect is 
more important for DSBs. In contrast, the actual fusion of 
the clumps (zc-zu) tends to make the n value of the bulges 
larger. 

Therefore we find that the physical encounter of the 
baryonic cores tends to produce bulges with larger shape pa- 
rameters (i.e. more concentrated), and n ~ 1 profiles tend to 
be formed by secular evolution. These results together with 
the fact that at z = most of GLOs show n < 1.5, regard- 
less of their past history, suggest that successive mergers 
have driven the bulges from exponential to more concen- 
trated profiles, to come back to exponential ones, if there is 
available gas in the systems. 

We think that the possibility of having secular evolution 
and core fusion associated to a merger helps to drive this cy- 
cle in the shape parameter. Note that the collisionless merger 
of stellar discs will not produce the same distributions. As 
a matter of fact, in order to detect this morphological loop, 
simulations should include gas dynamics and star formation 
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Figure 5. Disc scalelength (rj) as a function of bulge scalelength (r'b) for simulated objects at 2 = and progenitors (filled pentagons) 
for which bulges are described by < n < 5 (a), 0.5 < n < 2 (b), 0.8 < n < 1.5 (c) and 0.9 < n < 1.1 (d). Observational data from 
CdJB96 (crosses), MGH99 (open squares), KWKOO (open triangles) and MCH02 (open circles) for the same shape parameter intervals 
are also shown. 



since, otherwise, bulges can not be re-shaped back to n « 1 
(AOl). 

Fig. 1^ shows the variations of the n parameter during 
the ODP (a) and those produced during the fusion of the 
baryonic clumps (b) as a function of Mg^tarZ-^GTr*! where 
Mftar is the stellar mass formed in the progenitors at 2 and 
Mftar* is the total stellar mass of the GLOs at z = 0. This 
ratio is an estimate of the presence of a stellar bulge since in 
these simulations the stars form preferentially in the dense 
regions. T02 found that the properties of the potential well 
can be directly linked to the stability of a system, in the 
sense that the shallower the potential well, the more sus- 
ceptible the GLO to experience early gas inflows (i.e. sec- 
ular evolution). The formation of stellar bulges has been 
found to inhibit such early gas inflows. This picture is in 
agreement with our results, where the major changes in the 
shape parameter during a merger event correspond to the 
GLOs that experience DSBs {Mit^^/M^t^^ < 0.45), which 
are less stable objects and can be strongly perturbed by 
the collision with a satellite. In contrast, more stable sys- 
tems (Ms^tar/Afstar' > 0.45) would uot experience such im- 



portant perturbations in its mass concentration during the 
ODP (Table |). 

We can also see from Fig. ^a) that the most of the 
changes in n are negative and can be linked to smaller 
Mlt^JM^rJ values {M^t^JMICj < 0.45). This could be re- 
lated to the fact that stellar bulges can provide stability to 
the systems, making them less vulnerable to the influence 
of a satellite, principally during the ODP. From Fig. |^(b) 
we see that the fusion of the baryonic cores produce more 
positive changes in the n parameters and hence, a trend to 
increase it. Again, the larger variations are related to smaller 
bulges (i.e. smaller Mit^JMit^°). 

We also analysed the possibility that the changes in the 
n value during mergers are linked to the relative mass of the 
colliding systems, finding no correlation in contrast to the 
results of AOl. However, these authors only considered col- 
lisionless mergers while in all our mergers dissipation plays 
an important role. 

In Fig. 1^ we have plotted the averaged n parameters 
of the progenitors at different z intervals. We can see that, 
although the bootstrap errors are high, there is a trend for 
the averaged n to decrease toward n = 1 for decreasing z. 
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Figure 7. Variations in the shape parameter of the progenitors during the orbital decay phase (a) and during the fusion of the baryonic 
clumps (b) in single (filled circles) and double (filled squares) SBs, as a function of the fraction of stars already formed in the progenitors, 
^starZ-'^star ^hen the mergers start. M^^^^^ / = 0.45 (dashed lines) has been taken to be an indicator of the presence of a well-formed 
stellar bulge after T02. 




Nevertheless, at all redshifts we found the mean n smaller 
than n — 2, indicating that, in spite of the spread of the 
parameters, the majority of the bulges are best fitted with 
n< 1.5. 

3.3 The rb/rd ratio 

It has been argued by some authors (e.g., CdJB96) that 
there is a restricted range of r^/rd for spiral galaxies which 
has been interpreted as the result of secular evolution of 
the central mass concentrations. In Fig. ^ we showed the 
rd versus rb distribution for GLOs at 2: = and for 
observations. We have combined the observational results 



Figure 9. Ratio between bulge and disc scalelengths as a func- 
tion of the bulge shape parameter for simulated objects at z = 
(filled triangles), single SBs (filled circles) and double ones (filled 
squares). Observational data from KWKOO (open triangles) and 
MCH02 (open circles) are also shown. The minor box shows the 
changes of the plotted parameters for single (thin line) and dou- 
ble (thick line) SBs during the orbital decay phase (first arrow) 
and the fusion of the baryonic clumps (second arrow). 

from MGH99, KWKOO and MCH02 which include spiral 
galaxies from Sa to Sc. For these distributions we find 
< rb/rd >= 0.07 ± 0.02 (the boostrap error is aht = 0.02) 
and < rb/rd >= 0.10 ± 0.08 (dbt = 0.01) for the simulated 
and observed data, respectively, showing a good agreement. 
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However, we have already shown in Fig. ^ that the sim- 
ulated and observed scalelength distributions show a clear 
trend only for those objects with bulges with 0.9 < n < 1.1. 
Hence, we divided the observations and simulated data into 
three subsamples according to their n parameters {n < 0.9, 
0.9 < n < 1.1 and n > 1.1) and estimate their mean rb/rd. 
Table ^ summarizes the results. First, we note that obser- 
vational and simulated < rb/rd > in the three subsamples 
are in very good agreement, supporting the hypothesis that 
GLOs in CDM models have observational counterpart. We 
also see that the larger dispersion is found for bulges with 
n < 0.9 in both the observed and simulated data, and that 
these subsamples show higher < rb/rd > than the subsam- 
ples with n « 1, suggesting that these systems systemati- 
cally have larger rb (for this subsample: < rb 0.89±0.10 
and < rd >= 4.77 ±0.40). 

On the other hand, we find that systems with bulges 
better fitted by n > 1.1 show the smallest < rb/rd >. Hence 
for a given disc scalelength, the larger the concentration of 
the bulges, the smaller its scalelength (for this subsample: 

< rb >= 0.23 ± 0.03 and < rd >= 5.52 ± 0.30). If we as- 
sume that the n parameter correlates with morphology as it 
is indicated by several works (e.g., de Jong 1996; Graham & 
Prieto 1999, hereafter GP99; GdBOl), then the scalelength 
ratios should correlate with n. Fig. B shows the rb/rd ratio 
as a function of the shape parameter for the simulated ob- 
jects during mergers and for z = 0. We have also included 
observational data from KWKOO and MCH02. We find an 
anticorrelation between these two parameters, in agreement 
with observations. Assuming that the shape parameter is an 
indicator of Hubble type, this anticorrelation indicates that 
late-type spirals (n < 1.5) have larger rb/rd ratios than 
early-type spirals (n > 1.5). We also note that the relation 
for the observed and simulated ratios seems to flatten for 
n <C 1. A simple extrapolation gives a value of rb/rd « 0.40. 

Because the shape parameter and the scalelength of the 
Sersic law are coupled parameters, this trend of rb /rd with 
n could be the result of the fitting formula. However, Tru- 
jillo, Graham & Caon (2001, hereafter TGCOl) argued that 
there is a physical trend stronger than that resulting from 
the fitting formula. Because we use the same expression for 
all fittings, the fact that there is a segregation in the scale- 
length distribution according to the value of n suggests that 
this trend could have physical basis and could be related to 
a difference in the evolution of the GLOs or in the efficiency 
of some physical process. This is confirmed when r^ff is used 
instead of rb for which we find < reff/rd >= 0.29 (abt=0.04) 
for n < 0.9 , < rcft /rd >= 0.26 (abt=0.03) for 0.9 < n < 1.1, 
and < reff/rd 0.30 (crbt=0.05) for n > 1.1 (the com- 
bined observational sample of KWKOO and MCH02 shows 

< rcff/rd >= 0.23 (crbt=0.02) for n < 0.9, < r^ff/rd >= 0.28 
(crbt=0.04) for 0.9 < n < 1.1, and < r^s/rd >= 0.35 
(crbt=0.03) for n > 1.1). We have used the approximation of 
Ciotti & Bertin (1999) to transform simulated rb into reff. 

Let us now look at the mean < rb/rd > during merg- 
ers. In Table ^ we show the mean rb/rd at za and zo for 
SSBs and DSBs. We find no differences between the mean 
rb/rd but both samples agree in increasing the mean ratios 
and the dispersions at zu- Fig. ^ shows the distributions 
of scalelength during mergers. Note that, independently of 
the stage of evolution of the mergers, the distributions have 
similar characteristics and are within the observed range at 



Table 3. Mean values of the ratio between bulge and disc scale- 
lengths for simulations (SIM) and observations (OBS) for differ- 
ent cases. S and D denotes single and double SBs, respectively. 
The bootstrap errors are shown in brackets. 





< rb/rd >SIM 


< rb/rd >OBS 




O.lOiO.lO (0.01) 




2 = 


0.07±0.02 (0.02) 


0.10 ±0.080 (0.008) 


n < 0.9 


0.24 ±0.20 (0.05) 


0.16 ±0.07 (0.01) 


0.9 < n < 1.1 


0.13 ±0.06 (0.02) 


0.14 ±0.05 (0.02) 


n > 1.1 


0.05 ±0.05 (0.01) 


0.03 ±0.04 (0.01) 


S ZA 


0.07±0.05 (0.02) 




S 2D 


0.11±0.08 (0.03) 




D ZA 


0.07 ±0.08 (0.02) 




D 2D 


0.10±0.10 (0.04) 





z = 0. Finally in Fig. ^ in the small box, we show how the 
averaged r^^/vd moves throughout mergers associated with 
single (thin line) and double (thick line) starbursts. It is clear 
that early gas inflows during the ODP tend to decrease n in 
both cases but for DSBs the changes are more important, 
as expected from their stability properties (see Fig. ^ and 
T02). Conversely, the fusion of the cores moves the distri- 
bution toward higher < n >. 

3.4 The B/D ratio 

The luminosity B/D ratio has been traditionally used as an 
indicator of morphology (e.g., Andredakis & Sanders 1994). 
The correlation between B/D and n has been shown by 
APB95 and recently confirmed by TGCOl. Linear regres- 
sion through observations show d{B/D)/dn — 1.0 ± 0.3 
((Jbt = 0.3) for APB95, 0.8 ± 0.5 ((Jbt = 0.5) for KWKOO 
and 1.7 ±0.3 (abt = 0.3) for MCH02. The combined sample 
has d{B/D)/dn = 1.2 ±0.1 (abt = 0.1). 

Note that in the models we have masses instead of lu- 
minosities and that the mass-to-light ratios of the discs (7d) 
and bulges (7b) could be different because of their differ- 
ent stellar populations. However, their correct values and/or 
possible dependence on radius as well as redshift are still 
unclear. In order to match the luminosity B/D we have to 
rescale the mass B/D ratio by a constant factor of 20. We 
address that we are using the same mass-to-light ratios for 
all GLOs regardless of the redshiftj^ disc scalelength or ob- 
servational band. 

In Fig. ^ we show the B/D ratio for the simulated 
GLOs during merger events at the four redshifts of inter- 
est and the observational data from APB95, KWKOO and 
MCH02. The low statistical number restrict the use of lin- 
ear regression through the simulated data. However, we see 

^ Nevertheless, we found that very similar 7b /7d ratios are 
needed to match the observed luminosity B/D for progenitors 
at different 2, although our dispersion is quite high. Note that we 
have used the total baryonic masses to estimate luminosities. In 
a more realistic model including energy feedback, some fraction 
of the cold gas would be reheated. Hence, the values required for 
our models to match the observed luminosities could be taken as 
upper limits to the mass-to-light ratios. 
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Figure 10. Distribution of bulge and disc scalelenghts for the simulated galactic objects at za (a), zg (b), zq (c) and zo (d). Different 
symbols have been used to distinguish between n < 1.5 (encircled symbols) and n > 1.5, and between double (filled squares) and 
single (filled circles) starbursts. Observational data from MGH99 (n = 1 crossed squares, n = 2: open squares), KWKOO {n < 1.5: 
crossed-triangles, n > 1.5: open triangles) and MCH02 {n < 1.5: crossed-circles, n > 1.5: open circles) are also shown. 



that, at any time during merger events, the simulated values 
are within the observational range. 

In Figjl^ we plot the variations . — Wz^. versus B /D^^ — 
B/Dzj where i/j can be A, C or D. From this figure we see 
that GLOs which experience gas inflows show more impor- 
tant changes in both n and B/D during the ODP. In the 
case of those GLOs experiencing SSBs (e.g., no early gas in- 
flows) there is no important change in the shape parameter 
or the B/D ratios during the ODP. We also note that the 
changes in n and B/D do not correlate. Conversely, during 
the actual fusion of the objects (from zc to zo) both types 
of GLOs experience large changes in both parameters in the 
expected direction: the larger the increase in the shape pa- 
rameter, the larger the corresponding increase of the B/D 
ratio or vice- versa. Fig. |l^(c) shows the overall changes over 
the merger events. This distribution has higher dispersion 
than that obtained during the fusions because it also com- 
bines the effects of secular evolution. However, the overall ef- 
fects of mergers is to settle the correlation between B/D and 
n so that more concentrated bulges (higher shape parame- 
ters) have the larger B/D ratios. However, we note that, 
although the range of values can be matched, the simulated 



B/D do not show the same correlation signal (r = 0.20) 
than the observations (r = 0.55). 



4 CONCLUSIONS 

We have studied the properties of the baryon distributions 
in galaxy-like objects focusing on the effects of mergers. We 
resort to observations of spiral galaxies of different morphol- 
ogy to constrain our findings. 

We found that on average, galactic objects formed in 
hierarchical clustering scenarios reproduce the angular mo- 
mentum and structural parameter distributions of spiral 
galaxies, if a stellar bulge is allowed to form and early gas 
depletion is avoided. We have succeeded in these two aspects 
but at the expense of inhibiting star formation on discs. A 
consistent treatment of energy feedback may help to remove 
this caveat. 

These simulations have allowed us to study how mergers 
change the distribution of baryons by analysing the evolu- 
tion of their structural parameters. We found that, on aver- 
age, galactic objects tend to have nearly exponential bulges 
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Figure 11. Luminosity bulge-to-disc ratio (B/D) as a function of the shape parameter of the progenitors during single (filled circles) and 
double (filled squares) SBs at 2a (a), zb (b), zc (c) and zd (d). Observational data from APB95 (asterisks), KWKOO (open triangles) 
and MCH02 (open circles) are also shown. 



at all redshift. However, note that these simulations produce 
bulges with shape parameters in the range 0.5 — 4. 

For those systems with n « 1 bulges we found a corre- 
lation among their bulge and disc scalelengths in very good 
agreement with observations. However, for n > 1 the scale- 
length distribution is disorder and displaced to smaller bulge 
scalelenghts. Observations show the same behaviour. The 
opposite distribution is found for systems with n < 1 bulges 
which have larger rb/ra values and larger dispersions. We 
found the disc scalelengths to be approximately indepen- 
dent of n parameters, so that this correlation implies that 
more concentrated objects have smaller bulge scalelengths 
(or large effective radius). We also found a ma^ximum rb/ra 
of w 0.40 for n — » 0. Higher numerical resolution simulations 
are needed to study the formation of such low mass surface 
profiles. 

In order to test our results for low resolution problems 
in the determination of the structural parameters of the 
galaxy-like objects, we followed Steinmetz & Miiller (1994) 
and used the bootstrap technique to estimate the effects of 
low particle number statistics within the objects. We calcu- 
lated an accuracy better than 25 per cent for the scalelengths 
and shape parameters. 



We found that gas inflows during the orbital decay 

phase tend to produce important changes in the mass dis- 
tributions generating n ~ 1 profiles, while the fusions of the 
baryonic cores tend to increase the n parameter. As a con- 
sequence, a morphological loop can be driven by mergers 
which might be responsible of triggering secular evolution 
as well as of the baryonic core fusions. The triggering of 
secular evolution is found to be linked to the presence of a 
stellar bulge so that systems with well-formed stellar bulges 
do not experience early gas inflows. Hence, the pace of this 
morphological loop could be regulated by the properties of 
the galactic central potential wells which are also affected 
by the merger history of the objects (see also Tissera & 
Dommquez-Tenreiro 1998 and Tissera 2000). 

We found that the simulated mass bulge-to-disc ratios 
are within observed range. It is also noted that during the 
ODP larger changes are observed in the B/D ratios of those 
objects that experience gas inflows. The changes during this 
period are, however, quite disorder. It is at the fusion of 
the baryonic clumps that changes in the shape parameters 
are correlated with changes in the bulge-to-disc ratio. In 
our simulations, the actual fusions are responsible of sig- 
nificantly increasing the mass concentration at the centre. 
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Figure 12. Variations in the shape parameters of the progenitors during the orbital decay phase (a), the fusion of the baryonic elunips 
(b) and the whole merger process (c) in single (filled circles) and double (filled squares) SBs as a function of the corresponding changes 
in the bulge-to-disc ratio, B/D. 



Hence, wo found that the fusion of the baryonic cores could 
be the process that determine the observed correlation be- 
tween the luminosity B/D ratio and the shape parameter or 
morphological type. However, we found no dependence on 
the relative masses of the colliding objects. 

Overall, our results indicate that the morphological 
properties of galactic objects are the result of their merger 
histories within a hierarchical clustering scenario. Based on 
the good agreement found so far with observations we sup- 
port the hypothesis of mergers as the main morphological 
driver along the Hubble sequence. Consequently, the partic- 
ular and detailed history of substructure aggregation could 
be a key point in the determination of the astrophysical 
properties of galaxies. 
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